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Abstract 


Using the approximate endpoint property, we describe a technique for exist- 
ing solutions of the fractional q-differential inclusion with boundary value 
conditions on multifunctions. For this, we use an approximate endpoint 
result on multifunctions. Also, we give an example to elaborate on our 
results and to present the obtained results by fractional calculus. 
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1 Introduction 


Fractional calculus and q-calculus are some of the significant branches in 
mathematical analysis. The field of fractional calculus has countless applica- 
tions (for instance, consider [2, 10, 32]). Similarly, the subject of fractional 
differential equations ranges from the theoretical views of the existence and 
uniqueness of solutions to analytical methods (for more details, see [4, 5]). 
There has been intensive development in fractional differential equations and 
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inclusion (for example, see [9, 13, 38, 36, 35, 37, 29, 20, 34]). Also, great 
attention was devoted to the fractional differential inclusions (for more in- 
formation, see [4, 7, 12, 30, 40]). For finding more details about elementary 
notions and definitions of fractional differential equations one can study, for 
instance, [6, 26, 31]. 


In this article, motivated by [8, 33] and among these achievements, we try 
to stretch out the problem in a sense for the fractional q-differential inclusions 
with integral boundary conditions, in conformity with the definition of the 
fractional Caputo type q-derivative of order a and the fractional Riemann— 
Liouville type q-integral. We state the basic definitions and some properties 
of fractional q-differential inclusion from [18]. For this purpose, we consider 
and discuss the inclusion problem 


“Deu(t) € W(t, u(t), w/(t),w’(t)), (1) 


under conditions u(0) + u(p) + u(1) = Ak fo(s, u(s)) ds and 


“Dju(0) + “Dyu(p) + “Dzu(1) = ia fo(s, u(s)) ds, (2) 


“Df u(0) + “DP u(p) + “DP u(1) = i, fils, u(s)) ds, 
where a € (2,3],0<q¢p,6 <1, 7 € (1,2), fi: JxR—- R for i = 1,2,3, 
are continuous functions, W : J x R® — P.,(R) is a multifunction, and 
“De is the fractional Caputo type q-derivative for t € J = [0,1]. The set 
of all compact subsets of R is denoted by P.p(R). Also, we look into the 
existence of solutions for the fractional q-differential inclusion problem on 
the multifunction W : J x R"t! + P(R), 


“Deu(t) € W(t, u(t), “D7 u(t),..., Dz" u(t)) , (3) 


with conditions - 
u!(0) + ayu'(1) = 744 (DF u(p), (4) 
u(0) + agu(1) = oy, Tul), 


where a € (1,2],0<q¢,p,7%1 < 1,@—~% € [1, 00) for alll <i<n, 


here n > 1 andteé€ J = (0,1). 


As before, we remind some of the previous works briefly. In 1910, the sub- 
ject of q-difference equations was introduced by Jackson [22, 24, 23]. After 
that, at the beginning of the last century, studies on q-difference equations 
have been appeared in many works, especially in [1, 3, 15, 28, 39]. An excel- 
lent account in the study of fractional differential and q-differential equations 
can be found in [21, 25, 26]. 
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2 Preliminaries 


We recall from [14, 25] some basic definitions, notation, and results of q- 
fractional calculus, which are needed throughout this article. In fact, we 
consider the fractional g-calculus on the specific time scale T = R, where 
T,, = {O} U{s: 5 = soq”}, for a nonnegative integer n, so € R, and q € Jo. 
Let a € R. Define [a]q = (1—q*)/(1—q) [24]. The power function (y—z)@ with 
n € No is defined by (y—z) — ae (y—zq*), forn > 1 and (y— 2) =1, 
where y and z are real numbers and No := {O}UN [1]. Also, foro € Rand g 4 
0, we have (y— 2) = v7 Tro (y—z9")/(y—zq?**). If z = 0, then it is clear 
that y(?) = y? [11] (Algorithm 1). The q-Gamma function is given by P(o) = 
(1—q)°-) /(1—q)?-, where o € R\{0, —1, —2,...} [24] (Algorithm 2). Note 
that, [g(o + 1) = [o]g¥g(c). The q-derivative of function g is defined by 


and D,[g](0) = lim;_.0 Dg[g](7), which is shown in Algorithm 3 [1]. Further- 
more, the higher order q-derivative of a function g is defined by Dj|g](7) = 
D,[D?~"[g]](7), for n > 1, where D?[g](r) = g(r) [1]. The q-integral of a 
function g is defined on [0, b] by 


CoO 


tar) = f " g(é) dg€ = 7(1— 4) Ss aor), 


k=0 


for 0 < r < b, provided the series is absolutely converges [1]. If 7 in [0,7], 
then 


T lee) 
/ a(€) dé = Tglg|(T) — Tgla|(r) = 1 — a) Soa [Ta(T a") — ralra*)) , 


whenever the series exists. The operator I? is given by I?[g](r) = g(r) 
and I? (g|(7) = I,{1?-*[g]](r) for n > 1 and g € C((0,7)) [1]. It has been 
proved that D,[I,[g]](7) = g(r) and I,[D,[g]](7) = g(7) — g(0) whenever g is 
continuous at T = 0 [1]. The fractional Riemann—Liouville type q-integral of 
the function g on Jo = (0,1) for o > 0 is defined by I}[g](r) = g(r) and 


for t € Jo [17, 10](Algorithm 4). The Caputo fractional q-derivative of a 
function g is defined by 
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= mio [o ge) el-7-Y ple [al(€) dag, (6) 


where t € Jo and o > 0 [17]. It has been proved that I¥[I?(g]](r) = 
19t”[g|(r), and D¢ {IJ [g]](r) = g(r), where o,v > 0 [17]. Algorithm 4 shows 
pesudo-code II [g](7). 

We say a multifunction G : J > P.,(IR) is measurable whenever for each 
real number y, the function t > d(y, G(t)) is measurable [16]. The Pompeiu— 
Hausdorff metric Ha: 2* x 2* — [0,00) on a metric space (X, p) is defined 
by 


H,(A,B) = max { sup p(a,B), sup p(A,0) 


where p(A,b) = infaea p(a,b) [19]. The set of closed and bounded of X, 
and the set of closed subsets of X are denoted by CB(X) and C(X), re- 
spectively. In this case, (CB(X),H,) and (C(X),H,) are a metric space 
and a generalized metric space, respectively [27]. An element z € X is 
called an endpoint of multifunction W : X — 2* whenever Wz = {z} [8]. 
Also, the multifunction W has the approximate endpoint property whenever 
infrex SUPyewz P(z, y) = 0 [8]. A function 6 : R > R is called upper semi- 
continuous whenever lim sup,,_,5, 9(An) < @(A) for all sequence {An }n>1 with 
An > 2 [8]. 


Lemma 1. [8] Consider an upper semi-continuous function @ : [0,co) > 
(0,00) such that 6(t) < ¢ and lim infi+.0(t — O(t)) > 0, for all t > 0. Also, 
Assume that (X,) a complete metric space and that W: X > CB(X) isa 
multifunction such that Ha(W(zx),W() < @(p(a, y)), for all z,y ¢ X. Then 
W has a unique endpoint if and only if W has the approximate endpoint 
property. 


3 Main results 


Right away, we are ready to state and prove our main results. Foremost, we 
make the adjacent one. 


Lemma 2. Suppose that v € C(J,R), that a € (2,3), that 0 < B,q,p < 1, 
that y € (1,2), and that f; : J x R > R, for i = 0,1,2, are continuous 
functions. The unique solution of the fractional q-differential problem 


“Dp u(t) = v(t), (7) 


with conditions (2) is given by 
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a i 1 
+ al fo(s, u(s)) ds — A [Ie'v(1) + IS v(p)] 
0 
1 
+ a(t) f  fi(ssu(s))ds-+ aa(t) [Ej o(2) + 1g?) 
0 


+ (by + as() [ g2(s, u(s)) ds 
+ (bg + ag(t)) [Ig-7v(1) + I¢-7v(p)] : 


where 


_ 3tDg(2— 8) — (p + 1)Pg(2 — 8) 
Sipt-F 1) 
_ (p+1)P4(2 — 8) — 80q(2 — Yt 
3(p'* + 1) 
_ -60"" ’ +10 (8-0 
6(p!-8 + 1)(p?- a 
3(p'-8 + 1)T4(3 


| 
T 


Gps led : 7 
6(p?-* + 1)T4(8- Ta C= 
Bip? + I)pe-7 + Ly 
314(3—1)P,(3— ae Pee @) 
6(pt-8 + 1)(p2-7 + IT 4(3 — B)’ 
2(p + 1)(p?-* + IT, (3 — 7) 4(2 — 8) 
Gp? + 1)(p?-7 + 1) g(3 — 9) 
(p? + 1), (3 = 7)(p' + 10, (3 = B) 
bi F + ip? ref) ° 
—_ (p? oH IT, (3 = 1)(p\-P + 1)P,(3 ==) 
G(pt PF 4-1) Gees (3 = 8) 
2(p + 1)(p?-8 + INT, (3 a y)P (2 = B) 
6(p!—8 + 1)(p2-7 + 1)P4(3 — B) 


— 6) 


by = 


Proof. As we have known that the general solution of (7) is 


1 t 
u(t) = | (t — qs)'°-Yu(s) dgs + co + cit + cot?, (9) 


P,(a) Jo 


where c; (¢ = 0,1,2) are arbitrary real constants [26, 31, 38]. Thus, 
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ti-B Qc9t?2-F 

epeu(t) = 2-8 v(t) 4 = | ; 

i) TO Op) Tye) 
Qcot?-7 
SDV) SIS) ee, 
as Cee 


Hence, we get 


u(0) + u(p) + u(1) = 3c9 + (1+ p)er + (1 + p*)e2 + I v(1) + I? v(p) 
and 


ph +1 2(p2-8 +1) 


“Dyu(0) + Deu(p) + “Dgu(1) =e £,2=p) ° T4389) 
+I2-Fy(1) + 12-F u(y), 
Q-— 
“Dy u(0) + “Dyu(p) + “Dj r(1) = ante 
qd 


ele etl) ee Op) 


By employing the boundary conditions, we have 


3co+ (1 +p)a+(1 pico = f fo(s, u(s)) ds 
— Iyv(1) — Ij v(p), 


po +l Wp? F +1) _ 
Gyo Sp). POS) =| 


—12-Fo(1) -12-Fo(p), 


ap" +1) _ S,2(Ss S 
ogy =| fa(s, (s))d 


=I ull) a op). 


This is a linear system of equations of triangular shape, having co, c1, and c2 
as unknowns. By a back substitution, we obtain 
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al Pee \ me 5 [is e(1) +150) 


Ty (2 — ete 
3(pi-F + Enis 


Pen 
© 8(pl-8 +1) 


i J. fols, u(s)) ds + bg [I2-7v(1) + I2-7o(p)] , 


a= 7 T,(2 = 6) - [ fisuts)ya 


[Ue v(1) + ISP v(p)| 


(p'-8 +1) 
C2 ‘ a- a— 
= pF FD [19 By(1) 1, - 
(p? ? + DT, (3 — “4 
PPT GTA A Oued 


\ (p?-* + 1)T,(3 — he a a-y a-y 
(p 1-8 4+.1)(p 27 Hr, (3 — B) [3 v( )+Ty v(p)] , 


Q= pe [ hlssuena 


q(3 *) a-y ay 
ee v(1) + 13-7 v(p)] . 


At once, we replace co, c1, and cp in (9) and find the solution u(t) as we 
stated. 


Assume that ¥ = C?(J) endowed with the norm 


||ul] = sup |u(¢)| + sup |u’(t)| + sup |w'"(t)]. 
ted ted teJ 


Then (4, ||-|]) is a Banach space (see [16]). For u € 4’, we define the selection 
set Sy. by the set of all v € L'(J) such that v(t) € W(t, u(t), u/(t), w(t) 
for all t € J. For the study of problem (1) and (2), we shall consider the 
following conditions: 


(C1) The multifunction W : J x R? — P.,(R) be is an integrable and 
bounded such that W(-,71,22,73) : J — Pep(R) is measurable for 
all x; € R; 


(C2) The functions f; : J x R— R be are continuous and map 9: [0,00) > 
[(0,co) be a nondecreasing upper semi-continuous such that 


lim inf(¢ — 6(t)) >0 


and O(t) < ¢ for all t > 0; 
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(C3) There exist m,mo,m1, m2 € C(J, [0, 00)) such that 


Ha (Wit, T1, 22, ©3), Wit, as 9) 


3 
1 
< =. i-2, 
< rea rea ud @2La 2) 


and 


[fi(t, 2) — f(t, 2')| < (t)-(|2 — 2"), 


1 
Ai +An+A3°7 


for allt € J, x, 2’, x;, 2, € R, where 


Ml oo Mo]|loo QI MI] 60 
i= |p lo HHmolleo , _ 2llrll 


q(a@ + 1) 3 31(a + 1) 
5Dq(2— A) llmalloo , 10P4(2 — B)I|mIlo0 
3 3, (a — 8 +1) 


+ 10(21,(2 — 8) +T4(3 — B)) 
, (eal ee ills) 
30'4(3 — AT q(a—7+]) 
_ filmlloo , 2P4(2—)Ilmlloc 7 7 
a= [Bayt Teco. + Ode 6) +189) 
T4(3— 7) (IlmallooPg(a — 7 + 1) + 2llmllec) 
‘ ( T@-Arya—y+) )| 
rn -| lImlloo Eo) mlets=1 5) zine). 
eRe 1) Ty(a-7+1) 


(C4) Multifunction N : ¥ — 2* is given by 


N(u) ={hEX|FvE Sw : A(t)=w(d}, 


for each t € J, where by applying the notation in (8), we have 
w(t) = Tpu(t) + =i fo(s, u(s)) ds — 3 [I v(1) + I¢v(p)] 
1 
+ a(t) f  fi(ssu(s))as-+ a(t) [Ej o(4) + 13° (9) 
0 


+ (b1 + a3(t)) [ fals,u(s)) ds 


+ (b2 + aa(t)) [Ie—7v(1) + 12-7 v(p)| . 
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Theorem 1. The boundary value q-differential inclusion problem (1) and (2) 
has a solution, whenever the multifunction N : ¥ + P(#) has the approxi- 
mate endpoint property and conditions (C1)—(C4) are hold. 


Proof. We demonstrate that the multifunction N has an endpoint, which 
is a solution of the problem (1) and (2). Because the multivalued map 
t F W(t, u(t), u(t), u’(t)) is measurable and so has closed values, it has 
measurable selection and so Sy, is nonempty for all u € ¥. At present, we 
show that N(u) C ¥ is closed for u € Y. Let u € ¥ and let {an}n>1 bea 
sequence in N(u) with up, — x. For each n € N, choose up € Sy, such that 


In(t) = Tfun(t) + a fo(s, u(s)) ds — F [Ue on (1) + In (p)| 


+ ai(t) i fis, u(s)) ds + ae(t) [Te Fon (1) + ean (p)] 


1 
+ (1 +aa(t)) | fals,uls)) ds 
0 
+ (bz + aa(t)) [To Vem (1) + Ie~7un(p)| ‘ 
It is noteworthy that {vn }n>1 has a subsequence that converges to some v € 
L*(J), because W has compact values. Again, we denote this subsequence 


by {Un}n>1. It is easy to go over that vu € Sw» and x,(t) tends to x(t), 
where 


x(t) = Tpu(t) + al fo(s, u(s)) ds — ; [I'v (1) + I?v(p)| 


+ ar(t) : fils, u(s)) ds + ae(t) [1e-F v(1) + I? -Fu(p)| 


+ (b1 + a3(t) ) f talssuioya 


+ (b2 + aa(t)) [[g- 701) v(p)] , 


for each ¢t € J. This implies that « € N(u) and so N has closed values. Since 
W is a compact multivalued map, it is easy to check that N(u) is a bounded 
set for all u € ¥&. Now, we show that 

Ha(N(u), N(v)) < A(|ju — of]). 


Let u,v € & and let hy € N(v). Choose wy € Syw,y such that 


hy(t) = Tywi(t +3 f fos s))ds— 2 [ls wi (1) + Tfwi(p )] 


+ ‘of fi(s, u(s)) ds + a(t) [UP wi (1) = Ie? w(p)] 
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(bi + a3(t co) ff se fo(s, u(s)) ds 


(bp + ag(t)) [2-7 wy (1) + 18-7 w1(p)] , 


for almost all t € J. Put 
Wut) = W (t, u(t), u'(t), u(t), Woy = W(t, v(t), 0’(6), 0" (0). 
Since 


fe = 1 
i 
Ha (Wa, Wow) <i, ie 


x 8 (|u(t) — n(t)| + lu’) — v'(t)| + lw'"(t) — o (ED), 


for all t € J, there exists w € Walt) such that 


~__m(t) 
Ay + Ao + Ag 
x 6 (|u(t) — u(t)| + |u'(t) — v'(t)| 
+|u"(t) — 0""(t))) 


for all t € J. Consider the multivalued map G : J > P(R), which defines 
the set of all w € R such that w satisfies in (10). Since w, and 


|wi(t) —w| < 
(10) 


1 
= 0 / / A i : 
7) mO(|u v| + ju’ — o’| + Ju Di kewread 
are measurable, the multifunction 
GO) Wu), 0), 0"), 
is measurable. Choose w(t) € W (t, u(t), u’(t), u(t) such that 


|wi(t) — we(t)| 
1 
< ——_—__ mit 
= Cpa eee Pu, 


x p (lu(t) — v(t)| + fu’) — o' | + |e") — oI), 
for all ¢ € J. Now, Consider the element hz € N(w), which is defined by 
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ho(t) = Tpwe(t) + 3 J fo(s, u(s)) ds — 3 [Ie we(1) + IS we(p)| 


+ax(t) | fils, u(s)) ds + a(t) [12-Pw9(1) + 12-Pw9(p)] 


(b: +a3(t)) ji fo(s,u(s))d 


(bz + a4(t)) [I~ wo(1) + 1F-%we(p)] , 


for allt € J. Thus, 


|hi(t) — ho(t)| < Tp |wi(t) — we(d)| 
tah | fo(s, v(s)) — fo(s, u(s))| ds 


+3 = [ite }wi (1) — we(1) | ‘wi (p) — w2(p)| ds] 


ao 


ite of Ifi(s,v(s)) — fa(s,u(s))| ds 
+ |a2(t)| [T2-F wr (1) — we(1)| 
+19-? |wi(p) — wa(p)|] 


1 

bi + a3(t)| ii ifa(a,(#)) — fala, tu(e))| de 
+ |be + aa(t)| be oe we(1)| ds 
+12~7|w1 (p) — we(p)| ds] 

1 lM) oo F IMolloo 
Ss Ay ey EAS eae) me 3 

_2ilmllo  , 5Pg(2 = B)Ilmafloo 
' 8l,(a+1) - 3 
_ 100, (2 — B)|lmlloo 

3D g(a — 6 +1) 


+ 10(2P,(2 — 8) + T,(3 — B)) 


Ty (3 = 7) (lmallool'g(a = 7 + 1) + 2I|mll00) 
‘ ( 30'4(3 — B)Tg(a—y+1) )| 


Ay 
= leo, 


[ni (t) — Ag (t)| < Tp" wit) — wa(t)| 


TAD iB) er 
te (p!-8 +1) [ue ?\w: (1) — w2(1)| 
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+19? |w1(p ) = wa(p)|] 
+ last Lats. 0l6) — fo(s,u(s))|ds 


+ |ae(t)| [Eg 7 Jews (1) — oo 


+19~7—"|w1(p) — we(p) 


plete 2 B)||Mlloo 
: Ce den 


< arate 


+ (2q(2- B) +T'4(3 - B)) 
ie2 hme aa — 7 +1) + 2llmlloo) 
«(4 Ty(3— Bl g(e—y+D )| 


A(||u — oll), 


~ Ay + Ay + A3 


and 
|ny (t) — hg (t)| < ae t) — we(t)| 


<D ffals.v(9) ~ fo(s,u(s))|ds 


=“ Y+1 
038 — ‘ Yay, (1) — w 
(p2-7 +1) [Ue | i(1 ) 2(1)| 
HS wile) — w2(p)|] 


1 Il] o0 
S A+ ha + hy cess 
_Pa(38 = 7) (IlmealloP ala = 7+ 1) + 2llmlloo) 
Tgla= 41) 


3 
a ere e oe vl), 

where 
¥ (p?-? + IT, (3 — 7) q(2 — 8) 

(p'-? + 1)(p?-7 + 1), (3 — 8) 
T4(3 — y)(p'? + 1)P4(3 — Bt 
pr? + 1)(p?-7 + YP g(3 — 8)’ 

(p?-? + 1) (3 — TP q(2 — 8) 
(p18 + 1)(p?-7 + 1) (3 — B)P g(a — ) 

T'4(3 — y)P4(3 — B)(p' 8 +: 1t 

poe De te ks Bay). 


Hence, 
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||h1 — ha|] = sup|hi(é) — ha(t)| + sup [hy (t) — h)(t)| 
ted ted 


+ sup |hy(t) — hg (t)| 


lla — yl[)(Ar + Ao + As) = A([le — yl). 


Therefore, it is easy to get that 
Ha(N(u), N(v)) < A(\|u — vf), 


for all u,v € XY. On the other hand, the multifunction N has the approximate 
endpoint property. By using Lemma 1, there exists u* € ¥ such that N(u*) = 
{u*}. Thus, by employing Lemma 2, u* is a solution of problem (1) and (2). 


At present, we investigate the existence of a solution for the fractional q- 
differential inclusion problem with integral boundary value conditions 


“De u(t) € W (t, u(t), “D7 u(t),..., “Dz” u(é)) , 
u'(0) + ayu'(1) = D7j_, “Dy‘u(p), (11) 
u(0) + agu(1) = Yi, Tju(p), 


for the multifunction W : J x R"*+! — P(R), where t € J, a € (1, 2], n > 2, 
0<¢,p.7%<1,a-% > 1 for alll <i<n, and 


4 Examples and algorithms for the problems 


In this part, we give a complete computational technique for solving prob- 
lems (1) and (2), such that it covers all the problems and presents numerical 
examples solving perfect. To this aim, we consider a pseudo-code description 
of the method for the calculated g-Gamma function of order n in Algorithm 2 
(for more details, see the link https: //en.wikipedia. org/wiki/Q-gamma_ 
function). 

Table 1 shows that when q is constant, the g-Gamma function is an in- 
creasing function. Also, for smaller values of x, an approximate result is 
obtained with a fewer values of n. It has been shown by underlined rows. 
Table 2 shows that the g-Gamma function for values q near one is obtained 
with more values of n in comparison with other columns. They have been 
underlined in line 8 of the first column, line 17 of the second column, and 
line 29 of the third column of Table 2. Also, Table 3 is the same as Table 2, 
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but x values increase in Table 3. Similarly, the g-Gamma function for values 
q near to one is obtained with more values of n in comparison with other 
columns. Furthermore, we provide Algorithm 3 that calculates (D7 f)(2). 
Here, we give an example to illustrate our first main result, which applies to 
the different values q in Theorem 1. 


Example 1. Consider the fractional q-differential inclusion problem 


9 i 1 
Dé u(t) € 0 100 sin u(t) + Too °° u'(t) a 
. |u!"(t) | 
"100 \14 u(t) 7}? 
under the integral boundary conditions 
3 nee 
u(0) + u(7) +u(1) = | 30 © u(s) ds 
and 
2 2 5 2 1 es?-1 
“Dj u(0) + “DZ u(z) + “Dz u(1) = fo Ssp- cos u(s) ds, (13) 
“Dj u(0) + De u(2) + Dz u(1) = fp 2+! cosu(s) ds, 


where t € J = (0,1],a= 2 B= 2, y= 3, and p = 3 in (1) and (2). Consider 
the map W: J x R® > P(R) defined by 


Wit ) 0 Z si + : cos : les 
@1,%2,%3) = — sinz r24 : 
alias “100° 100” 100 \ T+ ial 
Also, f; : J x R—- R are define by 
7 t?—-1 2-41 
fo(t, x) = 99 C8" fi(t,xz) = : cosxz, fo(t,x) = — COS &, 


and N : C?(J) > 2°) is defined by 


N(u) = {h € C?(J)| Ju € Swe A(t) = w(t), 


for all t € J such that 


19 
2 


1 .s?-1 
s a(t) | = cos u(s) ds + aa(t) a} 
0 


va) +103) 


+ (b1 +a3(0)) [ 3000 cos u(s) ds 
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7 Zz 3 
+ (bs + aa(t)) [#o(1) +2 0(9)) | 


where 
a(t) = sea(Gt= al) 
: 3((3)3 +1) 
alt) = seals) = TalS)t 
: 3(2)3 +1) 
ag(t) = —S(6a)? + WT alg )PalS)t + 3(2)* + YPo(§)PalG)0 
6(($)8 +1)((9)3 + YT 4(5) 
3\4 4 4 4 a 3\2 2 (14) 
ee 6((4)3 + UP alg Pals) = SP a(g)Pals)(Ca)* + 1)t 
6((3)2 +1)(($)8 + D4 (4) 
p, — HUA) + Whe $ aC $) — (CH)? + WEa($)(H4 + 0G) 
6((3)3 + 1)((3)3 + IT, (4) 
Ae (3)? + DT 4(4)((3)3 + DEAS) — 8((3)% + YP G(4)P (4) 
6((3)3 + 1)((3)3 + IT 4() 
Put m(t) = 2, mo(t) = &, mi(t) = $", mo(t) = 284, and g(t) = ¢. 
Then, we have 
Ages rae _Mmolloc , 2llmlleo , Ta(F)llmilloo , 100 ¢(§)|Illoo 
aed ie ee ee ee 3 Brace) 
102 g(F) +P a(S) Ta(9) es eta) 2 La) 
| 30 4(¥)V4(33) 
IImlJoo  2VG(F)|I™Mllo0 
wee r4(35) 
) 204) +0 4)rG) Seeger) + UPL) 
l(3)0' (93) 
4 IImlloo, Pal 3) (lmalloolg(Z3) + 2llml0) 
| P43) T,(# 


In following data of Tables 4 and 5, it is easy to check that 


3 
1 
Ha(Wé(t, ur, U2, us), F'(E, v1, V2, U3)) < he heen |wz — Ug|) 


and that 
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1 


; — f; < : 
| f;(t, w) Filtv)| = ie -pheey 


(t)p(lu— vl), 


fort € J, j =0,1,2. Because sup,¢n0) ||u|| = 0, we have 


inf sup uv] =o. 
weC?(J) fp, | I 


Thus, N has the approximate endpoint property. At present, by applying 
Theorem 1, the system of fractional q-differential inclusions (12) and (13) has 
at least one solution. 


5 Conclusion 


The q-differential boundary equations and their applications represent a mat- 
ter of high interest in the area of fractional g-calculus and its applications in 
various areas of science and technology. q-differential boundary value prob- 
lems occur in the mathematical modeling of a variety of physical operations. 
The end of this article was to investigate a complicated case by utilizing 
an appropriate basic theory. In this manner, we proved the existence of a 
solution for familiar problems of q-differential equations under three bound- 
ary conditions (1)-(2) and (3)-(4) on a time scale and showed the perfect 
numerical effects for the problem, which confirm our results. 
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Supporting informations 


Table 1: Some numerical results for calculation of Tq(x) with q = 4 that is constant, 
v= 4.5, 8.4, 12.7 and n= 1,2,...,15 of Algorithm 2. 


v2=4.5 r=84 x= 12.7 n w=45 r=84 e127 


2.472950  11.909360 68.080769 9 2.340263 11.257158 64.351366 
2.383247 11.468397 65.559266 10 2.340250 =11.257095 64.351003 
2.354446 =11.326853 64.749894 11) 2.340245 =11.257074 64.350881 
2.344963  11.280255 64.483434 12) 2.340244 =11.257066 64.350841 
2.341815 11.264786 64.394980 13 2.340243 = =11.257064  64.350828 
2.340767 11.259636 64.365536 14 2.340243 = =11.257063 64.350823 
2.340418  11.257921 64.355725 15 = 2.340243 =11.257063 64.350822 
2.340301  11.257349  64.352456 


CONBDTKRWNH!!S 


154 Ranjbar and Samei 


Table 2: Some numerical results for calculation of [g(x) with q = 5 5, 2, a =5 and 
n=1,2,...,35 of Algorithm 2. 


i". 45 q= 3 mG =3 a2 


3.016535 6.291859 = 18.937427 18 2.853224 4.921884 8.476643 
2.906140 5.548726 14.154784 19 2.853224 4.921879 8.474597 
2.870699 5.222330 11.819974 20 2.853224 4.921877 8.473234 
2.859031 5.069033 10.537540 21 2.853224 4.921876 8.472325 
2.855157 4.994707 9.782069 22 2.853224 4.921876 8.471719 
2.853868 4.958107 9.317265 23 2.853224 4.921875 8.471315 
2.853438 4.939945 9.023265 24 2.853224 4.921875 8.471046 
2.853295 4.930899 8.833940 25 2.853224 4.921875 8.470866 
2.853247 4.926384 8.710584 26 2.853224 4.921875 8.470747 
10 = 2.853232 4.924129 8.629588 27 = 2.853224 4.921875 8.470667 
11 2.853226 4.923002 8.576133 28 2.853224 4.921875 8.470614 
12 2.853224 4.922438 8.540736 29 2.853224 4.921875 8.470578 
13° 2.853224 4.922157 8.517243 30 = 2.853224 4.921875 8.470555 
14 2.853224 4.922016 8.501627 31 2.853224 4.921875 8.470539 
15 2.853224 4.921945 8.491237 32 2.853224 4.921875 8.470529 
16 2.853224 4.921910 8.484320 33 2.853224 4.921875 8.470522 
17 2.853224 4.921893 8.479713 34 2.853224 4.921875 8.470517 


CANaoRWNYH))S 


Table 3: Some numerical results for calculation of T'g(x) with « = 8.4, q= * 5 2 and 
n=1,2,...,40 of Algorithm 2. 
q=3 q=3 g=3 0 q=3 q=3 q=3 


11.909360 63.618604 664.767669 21 11.257063 49.065390 260.033372 
11.468397 55.707508  474.800503 22 = 11.257063 49.065384  260.011354 
11.326853 © 52.245122 384.795341 23 11.257063 49.065381 259.996678 
11.280255 —50.621828 336.326796 24 11.257063 49.065380  259.986893 
11.264786 49.835472 308.146441 25  11.257063 49.065379  259.980371 
11.259636 49.448420  290.958806 26 11.257063 49.065379 259.976023 
11.257921 49.256401 280.150029 27 = 11.257063 49.065379 259.973124 
11.257349 49.160766 273.216364 28 11.257063 49.065378  259.971192 
11.257158 49.113041 268.710272 29 11.257063 49.065378 259.969903 
10 =11.257095 49.089202  265.756606 30 =11.257063 49.065378  259.969044 
11) 11.257074 = 49.077288  263.809514 31 11.257063 = 49.065378 259.968472 
12 11.257066 49.071333 262.521127 32 11.257063 49.065378 259.968090 
13 11.257064 49.068355 261.666471 33 11.257063 49.065378 259.967836 
14) 11.257063 49.066867 261.098587 34 11.257063 49.065378  259.967666 
15 11.257063 49.066123 260.720833 35 = 11.257063 49.065378 259.967553 
16 =11.257063 49.065751  260.469369 36 =11.257063 49.065378  259.967478 
17 =11.257063 49.065564  260.301890 37 =11.257063 49.065378  259.967427 
18 11.257063 49.065471  260.190310 38 11.257063 49.065378  259.967394 
19 =11.257063 49.065425 260.115957 39 =11.257063 49.065378  259.967371 
20 11.257063 49.065402 260.066402 40 11.257063 49.065378 259.967357 
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Table 4: Some numerical results of ['g(a) in Example 1 with different values of ¢ by 
Algorithm 2. 


nlo=7 a=] a=} a= O=9f G=% A= a=7% 
q=3 
1 | 1.1174 0.9592 0.9559 0.9673 1.1055 1.13815 1.2199 1.5327 
211.1311 0.9505 0.9448 0.9500 1.0747 1.0990 1.1824 1.4805 
3 | 1.13856 0.9476 0.9411 0.9444 1.0647 1.0886 1.1703 1.4638 
4]1.1371 0.9467 0.9400 0.9425 1.0615 1.0851 1.1664 1.4583 
5 | 1.1376 0.9464 0.9396 0.9419 1.0604 1.0840 1.1651 1.4564 
6 | 1.1877 0.9463 0.9394 0.9417 1.0600 1.08386 1.1646 1.4558 
7 | 1.1378 0.9462 0.9394 0.9417 1.0599 1.0835 1.1645 1.4556 
8 | 1.1378 0.9462 0.9394 0.9416 1.0599 1.0834 1.1644 1.4556 
9 | 1.1378 0.9462 0.9394 0.9416 1.0598 1.0834 1.1644 1.4555 
10 | 1.1378 0.9462 0.9394 0.9416 1.0598 1.0834 1.1644 1.4555 
q=3 
1 | 1.1069 0.9743 0.9772 1.0122 1.2620 1.3087 1.4715 2.1039 
2 | 1.1877 0.9526 0.9493 0.9662 1.1655 1.2049 1.3437 1.8906 
3 | 1.1522 0.9426 0.9364 0.9453 1.1221 1.1583 1.2865 1.7960 
4] 1.1593 0.9378 0.9302 0.9353 1.1015 1.13862 1.2594 1.7514 
5 | 1.1628 0.9355 0.9272 0.9303 1.0915 1.1254 1.2463 1.7297 
6 | 1.1645 0.9343 0.9257 0.9279 1.0865 1.1201 1.2398 1.7190 
7 | 1.1654 0.9337 0.9249 0.9267 1.0841 1.1175 1.2365 1.7137 
8 | 1.1658 0.9334 0.9245 0.9261 1.0828 1.1161 1.2349 1.7111 
9 | 1.1660 0.9333 0.9244 0.9258 1.0822 1.1155 1.2341 1.7098 
10 | 1.1662 0.9332 0.9243 0.9257 1.0819 1.1152 1.2337 1.7091 
q= 4 
1 | 0.9665 1.1206 1.1787 14118 2.6441 2.8906 3.8168 8.5184 
2 | 1.0284 1.0602 1.0963 1.2516 2.1063 2.2761 2.9085 6.0237 
3 | 1.0710 1.0218 1.0443 1.1539 1.8020 1.9312 2.4107 4.7288 
4]1.1018 0.9954 1.0091 1.0891 1.6109 1.7160 2.1053 3.9658 
5 | 1.1248 0.9766 0.9840 1.0438 1.4826 1.5723 1.9040 3.4780 
6 | 1.1421 0.9628 0.9657 1.0111 1.38927 1.4718 1.7646 3.1483 
7 | 1.1555 0.9523 0.9519 0.9868 1.3275 1.3992 1.6647 2.9162 
8 | 1.1659 0.9444 0.9414 0.9685 1.2793 1.3455 1.5913 2.7481 
9 | 1.1740 0.9383 0.9334 0.9545 1.2429 1.3051 1.53863 2.6235 
10 | 1.1803 0.9335 0.9271 0.9436 1.2151 1.2743 1.4945 2.5297 
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Table 5: Some numerical results of Ai, Ag, and A3 in Example 1 with different values 
of q. 


m] Av Ap AB 
aL 
1 | 3.3364 1.1763 0.4559 0.2013 
2 | 3.3955 1.1988 0.4592 0.1979 
3 | 3.4147 1.2061 0.4602 0.1968 
43.4219 1.2088 0.4606 0.1964 
5 | 3.4240 1.2096 0.4608 0.1963 
6 | 3.4244 1.2098 0.4608 0.1963 
7 | 3.4246 1.2099 0.4608 0.1963 
8 | 3.4251 1.2101 0.4608 0.1962 
9 | 3.4251 1.2101 0.4608 0.1962 
10 | 3.4251 1.2101 0.4608 0.1962 


2.9820 1.0480 0.4467 0.2234 
3.1379 1.1076 0.4552 0.2127 
3.2160 1.13875 0.4593 0.2078 
3.2553 1.1525 0.4613 0.2054 
3.2754 1.1601 0.4623 0.2042 
3.2852 1.1639 0.4628 0.2036 
3.2898 1.1656 0.4630 0.2033 
3.2923 1.1666 0.4631 0.2032 
3.2939 1.1671 0.4632 0.2031 
3.2943 1.1673 0.4632 0.2031 


SCO OANOTKR WN HEH 


— 


1.8371 0.6109 0.3881 0.3526 
2.0805 0.7071 0.4077 0.3130 
2.2800 0.7853 0.4216 0.2868 
2.4430 0.8488 0.4319 0.2686 
2.5751 0.9001 0.4395 0.2554 
2.6823 0.9415 0.4454 0.2457 
2.7686 0.9748 0.4499 0.2385 
2.8380 1.0016 0.4534 0.2329 
2.8943 1.0232 0.4562 0.2286 
2.9392 1.0404 0.4584 0.2253 
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Algorithm 1 The proposed method for calculated (a — 6) 


Input: a, b, a, n, q 
1is¢l 
2: if n =0 then 
3: pel 
4: else 

5 for k = 0 to n do 

6 s<sx(a—bxa*)/(a—b+ q***) 

7: end for 

8 pea*x*s 

9: end if 

Output: (a — b)( 


Algorithm 2 The proposed method for calculated P'g(«) 
Input: n, q € (0,1), « € R\{0, -1,2,...} 
1: p+ 1 
: for k =0 ton do 
pe plage) 
end for 
T(x) — p/(1— q)*~* 
Output: T,(z) 


go wh 


Algorithm 3 The proposed method for calculated (Dg f)(x) 
Input: ¢ € (0,1), f(a), « 
1: syms z 
: if =0 then 
g + lim((f(z) — f(q* z))/((. — 92), 2,0) 


2 
3 
4: else 
5 
6 


a (f(x) — f(q* x))/((1 — qa) 
Output: (D,f)(2) 


Algorithm 4 The proposed method for calculated (I¢ f)(x) 


Input: qg € (0,1), a, n, f(x), x 
1s¢0 
: for i=0 ton do 
pyre ong sere 
si stpfxq« f(x*q’) 
end for 
Lg & (2% * (1-4) #5)/(T4(2)) 
Output: (Ij f)(z) 


OP OE OND. 


